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Abstract. Let / bo a newform, as specified by its Hecke eigenvalues, on a 
Shimura curve X. We describe a method for evaluating /. The most interest- 
ing case is when X arises as a compact quotient of the hyperbolic plane, so 
that classical g-expansions are not available. The method takes the form of 
an explicit, rapidly-convergent formula that is well-suited for numerical com- 
putation. We apply it to the problem of computing modular paramctrizations 
of elliptic curves, and illustrate with some numerical examples. 



1. Introduction 

1.1. Computing modular parametrizations. Let E be an elliptic curve over 
a totally real number field F. By a modular parametrization of E, we mean a 
surjective morphism 

(1) Ja,cX B ^E 

of abelian varieties over F, where Xb is the Shimura curve attached to an order 
in a quaternion F-algebra B (see §§2.2 and 2.3). The images under (1) of certain 
special divisors on Xb are called Heegner points (see §4.1), and provide some of the 
most versatile examples of analytic constructions of solutions to algebraic equations 
(see e.g. [15, 14, 16, 9, 7, 39]). 

It is of interest to compute (numerically) the curves Xb, the maps (1), and the 
images thereunder of special divisors. Methods for doing so have been given in 
many, but not all, cases of interest. When Xb has very low genus (e.g., genus 0), 
one can work with an explicit algebraic description of Xb (see [8, 35]); in general, 
such a model becomes unmanageable and one must resort to analytic means (see 
[4, §3.2]). Methods for computing Heegner points using p-adic analysis (see [3, 11] 
and references therein) apply when B ramifies at some finite place, but there are 
cases of interest in which B ramifies at no finite place, such as when F ^ Q has odd 
degree and E is everywhere unramified; in that setting, Darmon-Rotger [4, §3.2] 
describe as "completely open" the problem of computing Heegner points. Even 
when p-adic methods apply, it is of independent interest, and described as an open 
problem by Greenberg [12, p. 5], to obtain a uniform archimedean-analytic method 
for computing Heegner points. 

We are motivated by the latter problem of computing Heegner points, and more 
generally modular parametrizations (1), using archimedean analysis. 

1.2. Evaluating modular forms. It is classical that the map (1) is given by a 
certain holomorphic differential dz on Xb, in the sense that the composition 

Pic°X B (C) S JacX B (C) -> E(C) S C/A 

The author was supported by NSF grant OISE-1064866 and partially supported by grant 
SNF-137488 during the completion of this paper. 
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of the Abel-Jacobi map for Xb with the inverse Weierstrass parametrization of E 
should agree, up to isogeny, with the linear extension of 



A fundamental result of Eichler-Shimura (see [30, §7.4], [39, §1.4.9], [2, §10.3]) 
characterizes the weight two newform fg (up to a normalizing scalar) by expressing 
its Hecke eigenvalues in terms of the arithmetic of E over finite fields. 

Thus the problem of computing modular parametrizations reduces to that of 
evaluating a newform of known Hecke eigenvalues on a Shimura curve. This problem 
is of independent interest. For example, near the end of his paper Shimura Curve 
Computations, Elkies [8, §5.5] writes: 

The reader will note that so far we have said nothing about com- 
puting with modular forms on Shimura curves. Not only is this an 
intriguing question in its own right, but solving it may also allow 
more efficient computation of Shimura curves and the natural maps 
between them, as happens in the classical modular setting. 

We distinguish three senses in which one can "compute" the space of holomorphic 
modular forms of given positive even integral weight on a Shimura curve Xb- 

(I) as an abstract vector space. 
(II) as an abstract Hecke module. 

(Ill) as a concrete Hecke module realized as a space of differentials on Xb- 

Solutions to problem (I) are given by classical dimension formulas [30, §2.6] and to 
problem (II) by some classical algorithms involving Brandt modules [26], modular 
symbols [31], and their generalizations (see [5, 6, 37, 13, 17] and references therein). 
In this paper, we address the passage from a known solution of problem (II) to one 
of problem (III). 

1.3. The role of Fourier expansions. If F = Q is the rational number field, 
B = M 2 (Q) is the matrix algebra, and (for simplicity) A Ma (Q) = Xq(N) is the 
usual compactification of Fo(iV)\H for some natural number N, then X M2 ^ has 
the cusp oo at which / = /m 2 (Q) admits a Fourier expansion 



for some complex coefficients a/(w), where a/(l) is nonzero and explicitly normal- 
ized. The ratios a/(n)/a/(l), which depend multiplicatively on n, are determined 
explicitly by the Hecke eigenvalues of /, and hence, via the Eichler-Shimura re- 
lations, by computable invariants of E. For example, if p is a prime not divid- 
ing N, then a/(p)/a/(l) = p + 1 — #_E(F p ). Thus one may compute the map 
Jac Xm 2 (q) — > E by summing the series 






obtained by integrating (2) term-by-term (see e.g. [7, 22]). 1 



We suppress here discussion of certain technicalities, such as the choice of cusp at which to 
expand /. 
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If B ^ M2(Q), then Xb has no cusps (equivalently, the group T has no parabolic 
elements), so fs does not admit an expansion of the type (2). The absence of such 
an expansion makes it more difficult to evaluate fs] in fact, we are not aware of 
any successful attempts to do so prior to the work described in this paper. 

1.4. Summary of results. We will describe 

(A) a general method for computing the values and line integrals of modular 
forms (and their derivatives) on Shimura curves Xb using only archimedean 
analysis (see §5), and 

(B) an explicit, "ready-to-use" formula, obtained by applying the method (A), 
in the simple but already interesting case that B is an indefinite quaternion 
algebra over Q and Xb is attached to an Eichler order of squarefree level 
(see §3). 

We apply the formula (B) to several numerical examples in §4, so that §§3, 4 
and 5 follow roughly a "theorem-application-proof" paradigm. Applications of the 
method (A) to other settings, such as those arising from Shimura curves attached to 
quaternion algebras over totally real fields, will be taken up in a future paper. The 
method itself is general, and extends to arbitrary automorphic quotients attached 
to unit groups of quaternion algebras over number fields. 

Explicit formula. To sidestep the notational preliminaries needed to state formula 
(B) precisely, we illustrate it here in the simplest nontrivial case. Let A(z) :— 
qJ| neN (l — q n ) 2i — X^nsN °A («■)<?" be the Ramanuj an/discriminant function on 
the upper half-plane H := {z = x + iy : y > 0}. It satisfies the functional equation 
A(— 1/z) = z 12 A(z) and spans the space of cusp forms of weight 12 on SL2(Z). A 
famous theorem of Deligne implies that |oa(?t-)| < n 11 / 2 r(n) with r(n) the number 
of positive divisors of n. 

One can evaluate A on the positive imaginary axis z = iy by summing the series 

(3) A(zy) = ^a A (n)e- 2 ^, 

n<EZ 

which converges rapidly provided that y is not too small. 2 A very special case of 
Theorem 3.1 is a new formula for A(iy), which, while visibly "worse" than (3), has 
the virtue of extending to modular forms on compact Shimura curves for which 
analogues of (3) are not available: 

Theorem 1.1. Let yi and y 2 be positive reals. Then 

(4) ( yi y 2 ) 6 A(i yi )A(iy 2 ) = £ a A (ad - bc)e 12i9 W 12 (27rr) , 

ad— 6c>0 

where r £ and 9 £ M/27rZ are the polar coordinates defined by 3 

re ie = ayjy 2 /y\ + d^y 1 /y 2 + i(b/ ^Jy\y 2 - Cyjy x y 2 ) £ C x 
and W\ 2 is the rapidly-convergent infinite sum of Bessel functions 
W u (x) = 2- 11 n 12 {nxK n (nx) - K 12 {nx)). 

2 The latter condition is not serious, as the functional equation A(i/y) = y 12 A(iy) allows one 
to assume that y > 1. 

■^For notational simplicity, we suppress the dependence of r and 6 on yi,y2,a,b,c,d. 
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Figure 1. The "first" modular form on a compact Shimura 
curve: 

Im(z) 4 / 2 /0) for / : H — !• C in the one-dimensional space 
S4(Tq(1)), as depicted by SAGE's complex_plot command and eval- 
uated via the method of this paper. Here Tg(l) = (J2t=i^ e i) n 
SL 2 (K), where d = [1,0; 0,1], e 2 = [l/y/2, l/y/2; 3/y/2, -1/V2], 
e 3 = [1/2,1/2; -3/2, 1/2] , and e 4 = [0, ^2; 3y/2, 0] . 

It is true, although not completely obvious, that the RHS of (4) converges ab- 
solutely, and in fact rapidly; see Remark 3.3. 

The general formula given by Theorem 3.1 is no less explicit in principle than 
that of Theorem 1.1. It follows that one can compute modular parametrizations 
by compact Shimura curves in terms of integrals of expressions essentially of the 
shape (4) (see §4 for numerical examples). 

In Theorem 3.9, we give formulas for the values of Shimura-Maass derivatives of 
newforms, which may be used to compute Taylor expansions with respect to suitable 
local parameters. Such formulas should provide an effective tool for computing 
equations of Shimura curves, but we do not pursue such applications here. 

The formula (4) is similar to those arising from standard methods for computing 
values of L-functions, so one interpretation of the results of this paper is that it 
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is no more difficult to evaluate automorphic forms on quaternion algebras than to 
evaluate L-functions of logarithmically-comparable conductor. 

General method. We turn to sketching the general method (A), of which a more 
detailed account will be found in §5. Recall, from §1.2, that the problem under 
consideration is to recover the values of a newform of given Hecke eigenvalues on a 
Shimura curve (see §2 for definitions). The abstract possibility of such recovery is 
the content of the multiplicity one theorem (henceforth abbreviated Ml) of Jacquet- 
Langlands (see §2.3), so one may view the problem under consideration as that of 
giving a computationally effective realization of Ml. From this perspective, our 
main tool suggests itself naturally, as we now explain. 

Jacquet-Langlands proved Ml by establishing a correspondence between au- 
tomorphic forms on B x and on GL2, identifying scaling classes of newforms on 
compact Shimura curves over Q with those on congruence covers of SL2(Z)\H. 
Their proof involved a comparison of instances of the Selberg trace formula, and 
was non-constructive. 

Shimizu [29] , generalizing earlier work of Eichler, realized the Jacquet-Langlands 
correspondence analytically as a theta correspondence between GO(B) and GL2. 
Since GL2 forms admit Fourier expansions, Shimizu's result reduces our task to 
making that theta correspondence computationally effective. 

A key new ingredient, which may be of independent computational and theoret- 
ical interest, is a technique for computing Petersson inner products on finite index 
subgroups of SLi2(Z) that does not require knowledge of an explicit fundamental 
domain (see §5.3). We have in mind the natural generalization to the setting of 
automorphic forms on higher-dimensional quotients attached to GL2 over a num- 
ber field, where explicit fundamental domains (and multidimensional oscillating 
integrals thereon) seem prohibitively complicated. 

1.5. Other approaches. A precursor to the method of this paper was imple- 
mented in Jan 2011, described in the notes [25], and presented at the Arizona 
Winter School in Mar 2011. A preprint of Voight and Willis [34] describes another 
technique, quite different from our own, for performing computations on compact 
Shimura curves. Inspired by methods of Stark and Hejhal for computing Fourier 
expansions of Maass forms on SL/2(Z), the authors compute the Taylor expansion 
of a modular form on a compact Shimura curve in a suitable local parameter by 
solving for the linear conditions imposed by the automorphy and Hecke relations. 
Both techniques seem worth developing: for example, while the approach described 
by Voight and Willis is sufficiently flexible to extend even to noncongruence or 
nonarithmetic Fuchsian groups, the method described in this paper seems strong 
for certain applications, such as computing modular paramctrizations of elliptic 
curves, in which there are already well-developed methods for obtaining the Hecke 
eigenvalues. We thank John Voight for several lively discussions of such matters. 

1.6. Acknowledgments. The problem considered in this paper, as well as our 
general strategy, was conceived over the course of several fruitful conversations 
with Kartik Prasanna; it is a pleasure to thank him for numerous helpful discussions 
that have contributed substantially to our own understanding. We thank Andrew 
Snowden for helpful discussions during a visit to the University of Michigan in 
Jan 2011, at which the first computations of the sort discussed in this paper were 
carried out. We thank Henri Darmon, Daniel Disegni, Jon Hanke, Eren Kiral 
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Mehmet, Victor Rotger, William Stein, Carlos de Vera, and John Voight for helpful 
discussions and encouragement during the Arizona Winter School 2011, of which 
we thank the organizers for their support. We thank Nahid Walji for his helpful 
feedback on an earlier draft of this paper. 

2. Background and notation 

2.1. Summary. The remaining subsections (§2.2 — §2.4) of §2 aim to make precise 
the summary below, and may be skimmed or referred to as necessary 

Groups. For each squarefree integer D possessing an even number of prime factors, 
and for each positive integer V coprime to D, we fix a Fuchsian group 

Tg(N) =RH SL 2 (M) < SL 2 (M) 

arising from a choice of Eichler order R of level N in a quaternion algebra B /q of 
discriminant D (see §2.2) together with a sequence of inclusions 

JicB4 M 2 (R) 

(see Definition 2.1). We take Tq(N) — Tq(N) for convenience. 

Newforms. For each T < SL 2 (K) as above and each positive even integer k, we 
define the set -^(L) of (holomorphic cuspidal) newforms on T (see §2.3). 

Jacquet-Langlands correspondents. Let V — Tq(N) and T = T n (DN). We re- 
call the Jacquet-Langlands correspondence as a natural bijection J r fe(r')/C x = 
J r fc(r)/C x between scaling classes of newforms (see §2.4). We say that a pair of 
newforms /' G J-fc(r'), / G -7"fc(r) corresponding thereunder are compatibly normal- 
ized if f(z) — X)n<EN o,f(n)e(nz) with a/(l) = 1 and the Petersson norms of / and 
/' with respect to standard hyperbolic measures coincide (see Definition 2.2). 

2.2. Quaternion algebras and Eichler orders. We collect here some back- 
ground on quaternion algebras and Eichler orders that will be relevant for what 
follows. We refer to [33] for details and proofs. 

Quaternion algebras over general fields. Let F be a field. A quaternion algebra B 
over F is a simple F-algebra with center F and dimension 2 2 = 4 over F. A basic 
example is the algebra Af 2 (F) of 2 x 2 matrices. A quaternion algebra is called split, 
or said to split, if it is isomorphic to A/ 2 (-F) (as an -F-algebra); it is called non-split, 
or said not to split, otherwise. For a, b £ F* , let B — (a, b\F) denote the F-algebra 



F © Fi © Fj © Fij, where i, j satisfy the relations i 



b, and ij 



-J*- 



Then (a, b\F) is a quaternion algebra over F, and every quaternion algebra over F 
is isomorphic to one of this form provided that F is not of characteristic 2. 

Let E be an extension field of F . The most relevant examples are quadratic field 
extensions of F , and completions of F if F is a number field. Say that E splits 
B, or that B splits at/over E, etc., if B ©^ E = M 2 (£?) splits. It is known that 
E splits B if and only if there exists an embedding B ^ Ajf 2 (i?) (of F-algebras). 
Every quaternion algebra over F is split by some quadratic field extension of F . 
For example, if a is not a square in F , then F(y^a) splits (a, b\F), and we have an 
embedding 



(a, b\F) 3 x a + x x i + x 2 j + x 3 ij \-> 



xo + xi^/a X2+xzyfa 
b(x 2 -x 3 ^/a) XQ-xi^/a 
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A quaternion algebra B comes equipped with a main involution l that induces 
the nontrivial automorphism of every quadratic field extension of F contained in 
B. If B M2(E) as above, then t may be obtained by restricting the adjoint map 
(cd) ^ (-c ~a)- We denote by 

det(a) = aa L £ F 
tr(a) = a + a'eF 

the (reduced) worm, and (reduced) trace on B. A The trace induces an F-valued 
bilinear form 

(h,b 2 ) =tx(hV 2 ) = b 1 b 2 + b{b 2 
on 5, called the trace pairing. 

Quaternion algebras over number fields. Suppose henceforth that F is a number 
field, and let \F\ be its set of places. If v £ |F|, say that v splits B, or that i? 
splits at/over v, if F„ splits B. For example, if a, & € F x , then (a, &|F) splits at v if 
and only if the Hilbert symbol (a, b) v , which is defined to be 1 if a is a norm from 
F v (Vb) and —1 otherwise, takes the value 1. 

Let denote the set all v £ |F| that do not split B. Class field theory implies 
the map 6 h> induces a bijection between the set of isomorphism classes of 
quaternion algebras over F and the set of finite even subsets of \F\; in the case 
B = (a, b\F), this amounts to the product formula lit/ ' = 1- 

The (reduced) discriminant Ub of B is the squarefree integral ideal composed 
of the primes at which B does not split. When B is a quaternion algebra over the 
rational number field F = Q, the discriminant Qb is an integral ideal in Z generated 
by a unique positive integer ds, which we shall also call the discriminant of B. 

Orders. An order R in B is a unital subring of B that is a lattice, i.e., for which 
dim^ R = diniQ B. A maximal order is an order that is maximal in B with respect 
to inclusion. For example, if o is an order in F(yfa) = F(BFi, i 2 = a, and b £ FHo, 
then o © oj is an order in (o, b\F); if o is an order in F, then M2(o) is an order in 
M 2 (F); if o is maximal in F, then M 2 {o) is maximal in M 2 (F). 

An Eichler order R C B is an intersection of two maximal orders in B. Let o 
be the maximal order in the number field F, and o v its closure in F v for each finite 
place v of F. The level of an Eichler order R C -B is the integral ideal A/" of F, 
coprimc to 0^, such that for each finite place v of F at which B splits, there exists 
an isomorphism B<g> F F v -4- M 2 (F V ) taking R<S> a o v to °^ ). When F = Q and 
N is generated by the positive integer N, we shall refer to as the level of R. 

Indefinite rational guaternion algebras. An indefinite rational guaternion algebra 
is a quaternion algebra over Q that splits at the unique real place oo of Q, or 
equivalently, for which there exists an embedding B M 2 {M) of Q-algebras. Class 
field theory implies that the map B i— > induces a bijection between the set of 
isomorphism classes of indefinite rational quaternion algebras and the set of positive 
squarefree integers having an even number of prime factors. 

Let B be an indefinite rational quaternion algebra. The map sending an Eichler 
order R C B to its level N £ N induces a bijection between the _B x -conjugacy 
classes of Eichler orders in B and the positive integers coprime to dg. For example, 



For each embedding B <-» M2(E), the reduced norm on B is the pullback of the determinant 
on M2(E). Thus the notation "det", while mildly nonstandard, should introduce no confusion. 
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if B = Ma(Q), then B is an indefinite rational quaternion algebra of discriminant 

dg = 1, and every Eichler order of (positive integral) level N in B is conjugate to 

I z z\ 
V iVZ z J ■ 

Definition 2.1. For each squarefree positive integer D having an even number of 
prime factors and each positive integer N coprime to D, we choose a subgroup 

(5) T°(N) <SL 2 (M), 
as follows. 

• If D = 1, we take T°(N) = T (N) = ( fa f ) n SL 2 (Z) = ( fa f) n SL 2 (K). 

• If Z? 7^ 1, we choose an indefinite rational quaternion algebra Z? of discrim- 
inant D, a real embedding B Ma(R), and an Eichler order i? C £> of 
level AT. We then let Tff(N) = R H SL 2 (K) be the image of the group of 
norm one units in R. 

The group Tq(N) is well-defined in the sense that its SL 2 (M)-conjugacy class is 
uniquely determined by D and N. The quotient Tq(N)\ SL 2 (R) is compact if and 
only if 1. 

2.3. Newforms on Shimura curves. We recall the notion of a (holomorphic) 
newform on a Shimura curve. 

Shimura curves. Let H = {x + : y > 0} be the upper half-plane. The group 
GL 2 (M) + of real 2x2 matrices with positive determinant acts on H by fractional 
linear transformations: ( " b d ) z := (az + b)/ (cz + d). 

Let T < SL 2 (M) be a lattice, that is to say a discrete subgroup of finite Haar- 
covolume. For example, Tq(N) is a lattice. More generally, one attaches a lattice 
r to each sequence of inclusions 

(6) R c B ^ M 2 (M), 

where B is a quaternion algebra with one split real place over a totally real number 
field, B M 2 (M) is a fixed real embedding (unique up to conjugation), and R C B 
is an order, by taking T — RD SL 2 (M) the be the image of the group of norm one 
units in R. 

For each lattice T, one defines a compact Riemann surface Xy by suitably com- 
pactifying the quotient r\H (see [30, §1.5]). By a Shimura curve, we mean an X-p 
for some V arising from a sequence (6). Such curves have canonical models over 
number fields (see §4.1 and [30, §9]). We call Xr a compact Shimura curve 6 if r\H 
is compact. 

Newforms. For each k € 2N, we recall the weight k slash operator : for each function 
/ : H — > C and each g £ GL 2 (M) + , the function f\/,g is given by 

dct(g) fc / 2 

flk9iz):= Jc^TW fi9) 9 = 

By an automorphic function of weight /c on T, we mean a smooth function / : EI — > C 
that satisfies f\k"f — f for all 7 G T. The space Sfc(r) of cusp forms of weight fc 
on r consists of those automorphic functions of weight k that are holomorphic and 
vanish at the cusps of T, if any. One knows that Sk(T) is finite-dimensional, and 
that ^(r) is isomorphic to the space of holomorphic 1-forms on Xp. 

^This definition is not intended to be exhaustive. 

^This terminology is a mild misnomer, because Xp is in all cases a compact Riemann surface. 
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Let r denote the set of all a E GL2(R) + for which TaT is a finite union of 
either left or right T-cosets. The Hecke algebra H(T) — C[r\r/r], which consists 
of formal finite C-linear combinations of double cosets Tar (a £ T), has a natural 
right action on Sfc(r) that linearly extends /|Tar = ^ f\k a j if TaT = |J Taj. By 
a newform in Sk(T), we mean an eigenfunction of T-L{T). 7 Let .Ffc(r) denote the set 
of newforms in Sk(X). The set .Ffc(r) is preserved under scaling by C x . 

Multiplicity one. For each newform / E J-fc(r), there is a character A/ : H(T) —> C x 
of the Hecke algebra with the property f\<p = A/ (<p)f for all <p E W(T). For general 
lattices T, a newform / need not be determined by A/; in other words, the one- 
dimensional irreducible constituents of the 'H(r)-module Sk(T) need not occur with 
multiplicity one. One perspective is that T is generally too small. 

For arithmetic lattices T arising from Eichler orders in quaternion algebras, one 
knows that T = B x n GL 2 (M) + is large, and (by a theorem of Jacquet-Langlands) 
that H(T) acts on Sk(T) with multiplicity one. 

In particular, each newform / on T = Tq(N) is determined by its character 
A/, the set J-k{T)/C x of scaling classes of newforms of given weight is finite, and 
the problem of recovering the values of a newform / from its Hecke data A/ is 
meaningful. 

2.4. Jacquet-Langlands correspondence. Let B be an indefinite rational quater- 
nion algebra, and let N be a positive integer coprime to d B . Set V = rg B (iV) 
and r = Tq(cIbN). The Jacquet-Langlands correspondence asserts (among other 
things) that there is a natural bijection of finite sets 



between the scaling classes of newforms on T' and those on T, characterized by 
a certain compatibility between the actions of H(T') and %(r). This map lifts 
non-canonically to a bijection ^(r') -H> ^.(T). 

Briefly, Jb € FkiT') corresponds to f E ^-(T) if and only if for each prime p 
not dividing dsN, one has A/ B (r'aT') = A/(raT) for some a' (resp. a) of integral 
trace and determinant p in V (resp. V). 

Definition 2.2 (Compatibly normalized Jacquet-Langlands correspondents). Let 
Jb € ^fc(r') and / E ^-(r) be Jacquet-Langlands correspondents. Write 



Call / and fs are compatibly normalized if a/(l) = 1 and the Petersson norms 
of / and fs with respect to the standard hyperbolic measures on r\H and r'\H 



coincide, i.e., J rXB y k \f(z)\ 2 ^ = J r ,^ y k \f B (z)\ 2 
Remark 2.3. It is known (see [28, eq. 22] or [33]) that vol(r\H) = | d B N l\ P \d B N( 1+ 
p- 1 ) and vol(r'\H) = fd B iV n p | dB (1-^)0^(1+ P" 1 )- 



Remark 2.4. Let / be as above, with a/(l) = 1. The coefficients a/(n) are deter- 
mined, in an explicit and straightforward manner, by A/, hence by A/ B . They are 
algebraic integers that may be computed in a variety of software packages (such as 



T k (T')/C x =T k (T)/C 



x 



(7) 




SAGE [32]). 



7 The concise definition given here, which is more closely aligned with the local theory of 
newforms after Casselman, agrees with the classical one when T = To(N). 



10 



PAUL D. NELSON 



Remark 2.5. By compatibly normalizing / and fg, we have determined the indi- 
vidual values only up to a scalar of magnitude one. However, for each pair 
of points zi, z 2 £ H, the quantity 

(8) ~hW)f B {z2) 

is well-defined in the sense that it is determined by / (hence by A/ , hence by A/ B ) 
and by the data R C B M 2 (R) giving rise to the definition of F'. Conversely, 
knowledge of the quantities (8) for all z%,Z2 £ H determines the form fg up to a 
scalar of magnitude one. 



3. Explicit formula 



In this section we state formulas for the values (see §3.1) and derivatives (see 
§3.4) of newforms on Shimura curves attached to squarefree level Eichler orders in 
indefinite rational quaternion algebras. §§3.2 and 3.3 introduce notation that may 
be referred to as necessary. 

3.1. Main formula for the values. Let B be an indefinite rational quaternion 
algebra, let N be a positive integer prime to cLb, and let R C B <^-» M 2 {M) be a 
sequence of inclusions giving rise to a lattice V = T^ B {N) as in §§2.2 and 2.3. Let 
L = Tq((1bN). Let k be a positive even integer, and let Jb £ .Ffc(r") and / £ ^(r) 
be compatibly normalized Jacquet-Langlands correspondents, as in Definition 2.2 
of §2.4. Let a/(n) be the nth Fourier coefficient of /, as in (7). Our immediate 
goal is to give an explicit formula for the well-defined quantities /b(zi)/b(^2) (see 
Remark 2.5) in the special case that N is squarefree. 

Thus, suppose also that N is squarefree. To state our formula, we must introduce 
some notation. Let d be a divisor of dsN . A basic consequence of Atkin-Lehner 
theory (see [1]) is that the Fourier coefficient a,f(d) does not vanish (in fact, that 
df(d) = ±d fe / 2_1 ), hence it makes sense to set 



(9) 



c f (d) 



n((d,d B )) fi(d) 



For each point z = x 



,1/2 





= /i((rf,AQ) 

d • a f{d) V d d • a,f(d) t 

iy £ H in the upper half-plane, define the matrix 



xy 



-1/2' 



V 



-1/2 



-1/2 



y 



-1/2 





'1 x~ 




y 




I 




l 



£ SL 2 (R), 



which has the property that it sends i to z under the usual action by fractional 
linear transformations. 

We postpone the introduction of some additional, more involved notation and 
R^ until §§3.2 and 3.3. Roughly, Vk(L,t) is a rapidly-convergent weighted sum 
over a lattice L < M 2 (R), while R^ is a lattice in B M> M2(K) that is locally 
dual to R at primes dividing d. Taking that notation for granted, our main formula 
reads: 

Theorem 3.1. For each Z\ = X\ + iyi and z 2 = x 2 + iy2 in H, we have 
(yiy 2 ) fc/2 7B^iT/ B (^) = ^a/(n) c f {d)V k (a^R^a Z2] n/d 



new 



d\d B N 



Proof. See Appendix A. 2. 



□ 
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Remark 3.2. Theorem 3.1 normalizes fg to have known Petersson norm (i.e., equal 
to that of /). This feature will be convenient in some expected applications. The 
Petersson norm of / may be computed in a variety of ways (see e.g. Example 5.7). 

3.2. Weighted sums over lattices. In this subsection we define Vk(L, t), which 
appeared in the statement of Theorem 3.1. Let k be a positive even integer, and 
let L be a lattice in M 2 (R). For example, M 2 (Z) is a lattice in M 2 (M). For each 
positive real x, define the infinite sum of Bessel functions 

W k (x) := 2 1 ~ k n k (nxK k -i(nx) - K k (nx)). 

n£N 

We refer the reader interested in how W k arises in nature to Appendix A. 2, and 
especially (25). 

There is a unique M-linear map i : C — > M 2 (R) sending i to ( _°i J ) ■ Let e = 
_ 1 J. Then e 2 = 1 and el(z)e = i{z) for all z G C. We have an orthogonal 
(with respect to the trace pairing) decomposition M 2 (R) = i(C) © ei(C), which 
gives rise to maps X, Y : M 2 (R) — > C by requiring that a — l(X(o)) + et(Y(a)) 
for all a G M 2 (R). As functions on M 2 (R), we have det = |A| 2 - \Y\ 2 . Set 
P := \X\ 2 + \Y\ 2 . Explicitly, 

vl , a + d + i(b-c) , . , a-rf + i(6 + c) , a 2 + fe 2 + c 2 + d 2 
-*(a) = 2 ' ^ = 2 ' ™ a ' = 2 

for a = G M 2 (R). 

Finally, for each positive real t, set 

(10) V k (L,t):= e lk ^ x ^W k (47r\X(a)\). 

det(a)=t 

Here e lka ^ x(a) = {X(a)/\X(a)\) k . 

Remark 3.3. To see that the sum defining V k (L,t) converges rapidly, note that W k 
is bounded and decays exponentially, and that for each x > 1, there are at most 
Ol(x a ) elements a G L for which det(a) > and |X(a)| < x. The latter claim 
follows from the inequality 2\X(a)\ 2 > 2|A(a)| 2 - det (a) = P(a). 

When computing (10) numerically, one restricts the sum over a to those elements 
of L for which the positive definite quadratic form P takes values smaller than some 
fixed large cutoff. 

Example 3.4. If L = M 2 (Z), then for each positive integer n we have 

V k (L,n)= Y e ike W k (2irr), where re 10 = a + d + i(b - c) G C x . 

ad—bc—n 

Remark 3.5. Unlike in the case of classical (parabolic) Fourier expansions, the spe- 
cial functions W k that arise in the definition (10), and hence in Theorem 3.1, are 
not representation-theoretically significant. Indeed, they are somewhat arbitrary; 
other choices are possible, although ours is convenient owing to its direct expression 
in terms of the classical if-Bessel function (see §5.3 and Appendix A. 2). They are 
instead analogous to the test functions that arise in the approximate functional 
equation method for computing values of L-functions outside their domain of ab- 
solute convergence (see [21, §5.2]). In fact, these is nonempty overlap between that 
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method and our own (see Example 5.7), although neither is a special case of the 
other. 

A generalization relevant for computing Shimura-Maass derivatives. The remain- 
der of this subsection will be needed only in the statement of Theorem 3.9. Define 
derivations 5±, 62 on the formal power series ring Z[[x, x, y, y]} via the table 



3 


S jX 








1 





-v 


—x 





2 





y 





X 



For nonnegative integers j\ and j 2 , define a polynomial Pj lt j 2 via 

Pj j — x~i 1+ i 2 e^ xS+v ^' '/ 2 S 3 ^ S 3 , 2 e~( xS+y ^ '/ 2 
The definition should be regarded as analogous to that of the Hermite polynomials 



H n (x) = (-l) n e* 



dx" 



An induction shows that Pj t j s has integral coefficients 



and is independent of x and x, so that we may write 

Pji ,h 



E 

ei,e 2 eZ> 



°ei,e 2 y V 



for some integers & e \ ,J e 2 2 . Now set 



(11) 



Vt J2 (L,t) 



E ^) {n+n+ei+e2)/2 C^^ e2 {L,t), 



where 

yh,h 

k;e 1 ,e 2 



.- e l{k+n+32)&rKX{oi) Y{a) ei Y{aY 2 W ] 



k+ 



31+J2 + e l+ e 2 



{Am\X{a)\). 



det(a)=t 



Remark 3.6. One can express Pj 1 J2 in terms of Laguerre polynomials 

(-xY 



(12) 



L^(x) = J2 



as Pji,hiy,y) = j-(-y) 32 3 y n ° L j V (yy) with, j = min(ii,j 2 ) and J = max(ji, j 2 )- 



l-i 



This identity is not essential for our purposes, so we omit the proof. 
Example 3.7. p . = 1, hence V°'°(L,t) = V k (L,t). p 0j2 = (-y) 12 . p jlt0 

y n - v\,h = (- 1 ) J2 y J2 ~ 1 (ro - h)- Ph,i = -y 3l ~ l {yy - h)- p 3 j - ; ' f ""' 



j\Lf\yy) 



t=yy 



3.3. Locally dual lattices. Let B be an indefinite rational quaternion algebra, 
and let R C B be an Eichler order of level N. We define here for each divisor d 
of dsN the auxiliary lattice R^ attached to R that appeared in the statement of 
Theorem 3.1 in §3.4. 

Let R* be the lattice dual to R with respect to the trace pairing. Thus, R is the 
set of all a G B for which (a, r) eZ for all r € R. We have R* D R. 

It is a general fact that lattices in vector spaces over Q may be characterized 
locally. In our context, this means more precisely the following: Let p be a prime 
number, and set B p = B ®z Z p = B <S>q Q p . For each lattice L < B (such as L = R 
or L = R*), set L p = L <E>z Z p . Then L p is lattice in B p . If L, L' are lattices in B 
that satisfy L p — L' p for all primes p, then L = V . 

We may characterize the dual lattice R* as follows: 
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• If p\d B N, then R p = R*. 

• If p | d,B, then B p is a discrete valuation ring, R p is the unique maximal 
order in B p , and R* is the inverse of its unique maximal ideal. 

• If p | N, there exists an isomorphism i : B p — > A'hiQp) taking R p to the 

Eichler order ( N \ £ ) . We then have ,(R* p ) = ( g ^ ) . 

For each divisor d of dgN, there is a unique intermediate lattice R* D R^ D R 
for which 

R p if p \ d, 
R* p if p I d 

as lattices in B p . We have RW = i? and = iT, while is indeed locally 

dual to R at the prime divisors of d. 

3.4. Extension to Shimura Maass derivatives. A generalization of Theorem 
3.1 gives a formula for arbitrary Shimura-Maass derivatives of a newform fg. These 
determine the geodesic polar expansion of Jb about an arbitrary point, which is 
useful for rapidly evaluating fg at a large number of points. 

Let k be a positive even integer, and T < SL 2 (IR) a lattice. Let 

be the Shimura-Maass raising operator sending (not necesssarily holomorphic) au- 
tomorphic functions on F of weight k to those of weight k + 2. This agrees with the 
action by the element | ( * J of the complexified Lie algebra of SL 2 (R) under the 
usual identification of automorphic functions / of weight k! (k! = k or k! = k + 2) 
on F with functions F(g) = tf\ k .g){i) on SL 2 (M) satisfying F (ig(™£ie SSl)) - 
F{g)e lk ' e for all 7 € T, 6> € K/2ttZ. For notational simplicity, we henceforth omit 
the subscript k from Sk and write 5 f for the £-fo\d composition Sk+2i-2°' ■ •°'5fc+2 '5fe- 

Example 3.8. If T contains ( 1 * ) and / € <Sfc(r) has the Fourier expansion 
f( z ) = SnGN a "'? n ' t nen its first couple of Shimura-Maass derivatives are Sf(z) — 
EneN (* - 47m) a„g" and S*f(z) = £„ eN " ^±f^ + (4™) 2 j a n q™ . 

For £ £ Z> , an induction confirms that S e f is the weight k + 2£ automorphic func- 
tion with Fourier expansion 5 f f(z) = y ~ e ^^2 ne f$ L^ k 1 \ATrny)a n q n , where ^ 
is a Laguerre polynomial as in (12). 

Theorem 3.9. With notation and assumptions as in Theorem 3.1, and Vi 1,j2 as 
in (11), we have 

y} +jl yj +j2 WM^j6^f B (z 2 ) = 5>,(n) £ c f {d)V^ (a^R^a^ 

neN d\d B N 

Proof. See Appendix A. 2. □ 

3.5. Complements. 

Specialization to CM points. The following brief remarks are not necessary for gen- 
eral comprehension of the paper. Let K be an imaginary quadratic field. Suppose 
that z = z\ = Z2 € H is the fixed point of some embedding K ^-s- B. Then the 
lattices a~ 1 R^a z arising in Theorems 3.1 and 3.9 decompose as "sums of trans- 
lates of products" of quadratic lattices tailored to the orthogonal decomposition 
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B = K K afforded by the trace pairing (see [27] for a more precise discussion of 
related issues), leading to a pleasant formula for <5 Ji /b(.z)<P 2 /b(;z) involving sums 
of the shape 

a /(^i TO i - hm2)r 1 (m 1 )r 2 (m 2 )W k | jl+ j 2+ei+e2 (4ir^/hT[). 

mi ,7712 >0 
fimi>£ 2 m 2 

Here £i and £2 are positive integers, while ri(m\) and r^im^) are the Fourier 
coefficients of certain imaginary quadratic theta series given by 

/ a \ k +Jl+J2 

rxirm) = X! (u|) ' r 2( m 2) = X! « ei « e2 

|a| 2 — mi I a| 2 =m2 

for some (translates of) lattices ai,a2 in K. We do not state any precise forms of 
such identities here, but remark that we have implemented their derivation algo- 
rithmicaUy. It would be interesting to compare what one gets in this way when 
ji = ji with what one would get via an explicit form of Waldspurger's formula and 
an approximate functional equation. 

In this way, one simplifies the problem of evaluating Vk(L; t) from one of enumer- 
ating short vectors in a rank four quadratic lattice over Q to one of enumerating 
short vectors in a pair of rank one Hermitian lattices over K. 

Remarks on assumptions. We have restricted N to be squarefree in our formulas 
only because our method of choice requires knowledge of the Fourier expansion of 
/ at each cusp of T (d B N) (see §5). The problem of writing down such expansions 
for forms of non-squarefree level, while not altogether trivial, seems orthogonal to 
the primary purposes of this paper (see e.g. [24, §1.9] for further discussion). For 
weights k > 2, one can use the holomorphic projection method described in §5.4 to 
obtain formulas valid for general levels involving the Fourier expansion of / only 
at the cusp 00. See Remark 5.8 for some musings on whether this method can be 
extended to the crucial case k = 2. 

4. Numerical examples 

In [25], we gave several examples of computations of the values and derivatives 
of modular forms on compact Shimura curves using a precursor of the method 
presented thus far. For example, we checked the rationality of their absolute values 
at CM points up to powers of known real periods. 

Shortly after [25] was presented, we improved the method to its present form and 
worked out some numerical examples along the lines of the motivation discussed in 
§1.1. In this section, we record several such examples. 

In more detail, this section is organized as follows. In §4.1, we recall how given a 
(modular) elliptic curve E/q and an imaginary quadratic discriminant D satisfying 
certain compatibility conditions, together with an auxiliary Atkin-Lehner operator 
t, one can attach a point P(D,r) G E(Hrj) (with Ho the ring class field attached 
to D) which is well-defined up to automorphisms of E and Galois conjugacy. We 
formulate a precise sense in which it is meaningful to compute such points (Problem 
4.1). In §4.2, we recall how such computation reduces to the evaluation of modular 
forms on Shimura curves. In §4.3, we record a small sample of numerical results 
obtained using the formulas stated in §3. 
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4.1. Setting. Let B be an indefinite rational quaternion algebra, let N be a posi- 
tive integer coprime to d B , and let R C B <^-> M 2 (IR) be an Eichler order and fixed 
real embedding giving rise to T' = T^ B (N) = Rn SL 2 (M) as in §§2.2 and 2.3. Let 
X B := X$ B (N) be the Shimur a curve attached to L' as in §2.3. 

CM points. Let D be a negative quadratic discriminant, so that D — oIkc 2 for 
some c £ N and some imaginary quadratic field K of discriminant < 0. 
There is a unique (up to isomorphism) quadratic order O — O b of discriminant 



D, constructed as O = Z[c(o!k + v / ^i<")/2]. Let Ho/K be the ring class exten- 
sion of K attached to O, so that class field theory gives a natural isomorphism 
Art : Pic(O) -> G&l(H D /K). Let h(D) = #Pic(C) = [H D : K] be the class 
number of D. 

Call an embedding O <—> R optimal if its extension to K ^ B satisfies KDR = O. 
Let CMo(ls) denote the set of images in X B of fixed points in H for optimal 
embeddings O i? (recall that i? M2(R) is fixed). There is a faithful action 
of Pic(C) on the set of optimal embeddings O R, which descends to a faithful 
action on CM D (X B ) that we denote by ( ■ a for C € CM D (X B ) and a € Pic(O): if 
C is fixed by an optimal embedding j : O <-» i?, and 7 G S n GL2(M) + is chosen to 
satisfy -fj(a)R — R, then ^ • a= 7C- 

Atkin-Lehner operators. Let 7? + = Bfl GL2(M) + be the monoid of positive norm 
elements in R. Call each r € i?+ that normalizes T' an Atkin-Lehner operator. 
The image in Aut(r'\H) of the Atkin-Lehner operators is a group that we denote 
by AL(r'). We have AL(r') = (Z/2) ul( - dBN \ where uj(d B N) is the number of prime 
factors of d B N, with representatives = T p indexed by the positive squarefree 
divisors d of d B N, where: 

• For p I d Bl t p e i? + is an arbitrary element of norm p. 

• For p a \ \ N,r p £ R+ is an element of norm p a that maps to ( «* 'q 1 ) under 



Canonical models. We regard Xb as an algebraic curve over Q with respect to its 
canonical model as in [30, §9], characterized as follows: for each negative quadratic 
discriminant D and each £ £ CM B (X B ), one has £ £ X B (Hfj); moreover, for each 
a £ Pic(0£>), one has C' a = £ Art ( a ). The group AL(r') acts on X B via maps defined 
over Q, hence preserving CM£>(X B ) and commuting with Gal(H B / K). The group 
AL(r') x G&\(H D /K) acts transitively on CM D (X B ). 

Heegner points. Let E be an elliptic curve over Q of conductor d B N, and let O = 
Od C K = Q(v^D) be an imaginary quadratic order of discriminant D. We assume 
that N is squarefree, that c is prime to N, and that CM B (X B ) 7^ 0. In view of the 
former assumptions, the latter holds if and only if 

(1) the prime divisors of d B do not split K, 

(2) the prime divisors of N split K, and 

(3) c is prime to d B . 

Note that these conditions determine d B in terms of E and K. 

The modularity theorem, the Jacquet-Langlands correspondence and Faltings's 
isogeny theorem imply that there is a surjective Q-map ip : JacXg — > E arising, 
as in §1.2, from some newform f B £ ^(r'). Such maps are not unique, since 



a surjection R — > A^Z/p* 3 ) taking R to 




for some /3 > a. 
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post-composing one with a Q-isogeny E — > E gives another, but we may and shall 
pin down the Aut(£')-orbit of ip by requiring that it not factor through a Q-isogeny 
E — > E of degree > 1. 

To a CM point £ £ GM-d{Xb) and an Atkin-Lehner operator r G AL(r'), we at- 
tach a point P(D, r) € E(H D ) obtained as the image of [rC] — [C] G Pic°(A B )(7J £) ) = 
Jac(As)(iJi)) under Although P(P/, r) depends upon (, its orbit under Aut(E) x 
Gal(PD/-Rr) depends only upon D and r. To classify such orbits conveniently, we 
assume that E is given by an affine equation y 2 + a±xy + 03 y = x 3, + 02a; 2 + 042; + olq 
with cocfhcicnts in Q (so that the {±l}-orbits on E are classified by x-coordinates), 
and also that E admits complex multiplication by neither Q(i) nor Q(v— 3) (so that 
Aut(.E) = {±1}). 8 The Galois orbit of an algebraic number is classified by its mini- 
mal polynomial. Letting x(P(D,t)) € Hp denote the ^-coordinate of P(D,t), the 
following computational problem is then meaningful: 

Problem 4.1. For (E, D, r) as above, find the minimal polynomial of x(P(D, r)). 

The novelty here is that we are able to address this problem via archimedean 
analysis even when 4^1 (see §1.1). 

4.2. Method. Our method for addressing Problem 4.1 is the natural one suggested 
by the content of §1.2 and Theorem 3.1. Let the newform f B € J-^F') realize 
ip : JacAs — > E, as in §1.2. Suppose that f B is compatibly normalized with its 
Jacquet-Langlands lift / as in Definition 2.2, so that Theorem 3.1 applies. The 
Fourier coefficients a f (n) may be read off from the arithmetic of E over finite fields 
in the usual way. 

(1) We find a basis for Hi(X B ,Z), which amounts to finding generators for (the 
abelianization of) V. In the examples below, we either found such genera- 
tors in the paper [23] or worked them out "by hand" . In more complicated 
examples, we could use the algorithm [36] as implemented in MAGMA. 

(2) Using Theorem 3.1, we compute to high precision some generators 9 for the 

lattice A fB := {/ 7 f B (z) dz : 7 G #i(A B ,Z)} = { f£° f B (z) dz : 7 G T'} 
(z e H fixed). 

(3) Using the built-in functionality of SAGE/PAPJ, we choose a Weierstrass 
parametrization C/Ae = E(C) for E. 

(4) We compute a minimal isogeny C/Af B —¥ C/Ae, or equivalently, a nonzero 
complex number fi for which (iAf B is contained in A# with minimal index. 
Our assumption Aut(P) = {±1} implies that /i is determined up to ±1. 

(5) For t G AL(T') and 71,72 € V, we compute f 1 ^ /b(^) dz via Theorem 

3.1 and its image P(C,r) under C -> C/A /s A C/A E = E(C) via the 
elliptic-exponential functionality of SAGE/PARI. Note that P(£,r) is in- 
dependent of 71 and 72, which we choose to minimize the length of the 
integration contour. 

(6) Let C u ...,Ch(D) be the Gal(P D /A)-orbit of some C € CM D (X B ). We 
compute these as the fixed points for the Pic(C)-orbit of some optimal 
embedding O R. The minimal polynomial of x(P(D,t)) divides J][(i — 
x(P(d,T))) e Q[t], whose coefficients we compute precisely. 



'The latter condition is imposed strictly for convenience (see [30, §6.1]). 
'in practice, it suffices to compute only one generator to high precision. 
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4.3. Results. Wc consider the following (Cremona-labeled) elliptic curves: 



E 


model 


15A 1 


y 2 + xy + y = x 3 + x 2 - lOx - 10 


21A 2 


y 2 +xy = x 3 - 49x - 136 




y 2 + y = x 3 + x 2 + 9x + 1 


26A ± 


y 2 + xy + y = x 3 — 5x — 8 


26B 2 


y 2 + xy + y = x 3 - x 2 - 213x - 1257 



The table below records some Heegner points that we found on these curves, all 
coming from uniformizations by compact Shimura curves {ds ^ 1)- The points not 
labeled as torsion have infinite order. 



E 


D 


h 


T 


minimal polynomial of x(P(D , t)) 


ID/il 


7 


i 


T 3 


x = — 4 








T~5 


ry> -< f ) f AVOl An 1 

X — O ^Z-LOIblUIlJ 








T~15 


X — ( 




7 . 9 2 
— / • Z 


i 
i 


7"3 


rr 1108/9^ c — 292977I 

X — — llUo/ZO — — Z Z 1 ( 




An 


o 
z 


T 3 


™ on /qo 

x — — zy / oz 








T~5 


X — O ^Z-LUISlUIlj 




A'X 
— *±o 


1 
1 


T 3 


x — — looy/zoo — — z 10 iuo 








7"5 


x = 3 (2-torsion) 








T 15 


™ 7Q/9107 7 — 2 A^i~ 1 7Q 1 




— oz 


z 


T 3 


ol0c2 2 (^7079^ J- 9091 77 
Z OX — Oiul ZX "T ZUZ1 / / 




-55 


4 


^3 


1 la; 2 + 2509a; - 751 




-67 


1 


7"3 


x = -461899/30976 = -2-*ll- 2 127 1 3637 1 


2\A 2 


-15 


2 


T7 


x = -23/5 




-39 


4 


T7 


13a; 2 - 103a; + 433 




-43 


1 


77 


x = -1588/225 = -3- 2 5- 2 2 2 397 1 




-51 


2 


77 


x = -116/17= -17~ 1 2 2 29 1 




-15- 2 2 


2 


77 


x = -527/5 = -5- L 17 L 31 L 




-67 


1 


77 


x = -17972/4225 = -5" 2 13- 2 2 2 4493 


35Ai 


-7 


1 


7-5 


x = -1/7 




-8 


1 


7-5 


x = -49/8 




-23 


3 


7-5 


23x a + 5062x 2 + 3951x + 41291 




-7-2 2 


1 


7-5 


x = -7401/7 = -7" i 3 i 2467 i 




-8-2 2 


2 


7-5 


x = -49/8 




-43 


1 


T 5 


x = -489/688 = -2- 4 43- 1 3 1 163 1 




-7-3 2 


4 


7-5 


7x 2 - 13x + 73 




-67 


1 


T 5 


x = -202489/4288 = -2- 6 67~ 1 7 1 28927 1 


26Ai 


-59 


3 


T 2 


ll 4 59 1 a; 3 + 2708827x 2 + 2483001x + 1553249x 




-67 


1 


7"2 


-2020489/2815675 = -5- 2 41- 2 67~ 1 2020489 1 


26B 2 


-24 


1 


Tl3 


-359/12 = -2- 3 3- 1 359 1 




-52 


2 


Tl3 


point at infinity 



We emphasize that although we have listed only the (minimal polynomial of the) 
x-coordinate, the y-coordinate may be recovered by solving a quadratic equation. 



Remark 4.2. These points were evaluated numerically using a mixture of SAGE/PARI 
and C++/GSL, and then recognized via continued fractions with the naked eye. 
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We do not claim that our results are correct, but it is likely that they are, be- 
cause in all cases we obtained a point with the expected field of definition. In 
principle, one could carry out such computations provably correctly by bounding 
the error terms in our explicit formulas and the heights of the points P(D,t) (via 
Gross-Zagier/Zhang formulas, see [22, §5]). 

Example 4.3. For (E,D) = (35Ai,-23) (resp. (26At,-59)), we have h(D) = 3, 
so there are three Atkin-Lehner orbits of CM points of conductor D on the Shimura 
curve Xq b (1) with ds = 35 (resp. 26). Letting xi,X2,%3 denote the x-coordinates 
of the resulting points on E(C), we compute numerically that 23 1 J\(t — Xj) (resp. 
11 4 59 1 Yl(t — Xi)) is approximately 



ll 4 59 1 t 3 + 2708826.9996 . . . t 2 + 2483000.9994 ...t+ 1553248.99991 . . . .) 

In each case we obtain approximately the integral-coefficient polynomial given 
above, whose roots generate a cubic subficld of the Hilbert class field Hd of Q(yl}). 

Example 4.4. We have listed only one example in which our computation returned 
a torsion point for "non-obvious" reasons: E = 26£>2, D = —52 and r = 7*13, for 
which the rank of the quadratic twist Ed is 0. 

Example 4.5. For E = 35A, D = —63 = —7 • 3 2 , we got a point defined over the 
ring class field H_ e3 = Q(\/=7, v 73 ^) of the order Z[3(l + \/-7)/2] C Q(\^7) of 
conductor 3. 



Throughout this section, we assume the notation and definitions of §2. 

5.1. Overview. Let B be an indefinite rational quaternion algebra, let N be a 
positive integer coprime to d B , let T' = T^ B (N), let R C B ^ M 2 (K) be the 
Eichler order and real embedding for which V = Rf] SL 2 (K) arises as (the image 
of) the group of norm one units in R, and let V — T^dsN). Let A: be a positive even 
integer, and let fg e ^fc(r') and / e ^fe(r) be compatibly normalized Jacquet- 
Langlands correspondents. 

One can compute the values of Jb in the following way: 

(I) Write down an explicit form of the Shimizu correspondence. This expresses 
fs in the form 



where Zl ,z 2 is an explicit non-holomorphic theta series depending only 
upon R C B M- M 2 (M) and k. Thanks to a result of Watson [38, §2.1], 
doing so amounts to an exercise in "unadelization." We define Zl ,z 2 m §5-2 
and show that it satisfies (13) in Appendix B. 
(II) Compute the Petersson inner product on the RHS of (13). We discuss sev- 
eral approaches to this problem in §5.4. Our preferred approach (Theorem 
5.6), which leads to the explicit formulas of §3, requires also the following 
intermediate step: 



23t 3 + 5061.9999982 . . . t 2 + 3950.9999999 ...t + 41290.9999609 . . . 



(resp. 



5. General method 



(13) 
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(l|) Write down the constant term in the Fourier expansion of y k f(z)0 Zl ,z 2 i z ) 
at each cusp of T. We carry this out in Appendix A.l when N is 
squarefree. 

After stating a precise form of the identity (13) in §5.2, we turn to a discussion 
of the general problem of computing integrals of automorphic functions (e.g., Pe- 
tersson inner products). This discussion occupies the remainder of §5 and serves to 
motivate Theorem 5.6, which gives a formula for rather general such integrals via 
roughly a "hybrid Rankin-Selberg and approximate functional equation" approach. 
In Appendix A. 2, we apply the procedure outlined above to prove the formulas 
stated in §3. 

Remark 5.1. This procedure applies also to more general automorphic quotients 
(such as Shimura curves associated to quaternion algebras over totally real fields), 
but we reserve further comments in that direction for a future paper. 

5.2. Definition of Shimizu theta function. Let notation be as in §5.1. Recall 
the functions X and P defined in §3.2. For each lattice L < Af 2 (IR) and positive 
reals t and y, set 

U k {L-t lV ) :=\y £ X{a) k e- 2 ^ p ^. 

det(a)=t 

Fix zi, z 2 G H, recall the notation a Zi from §3.1, and define 10 
9z u z 2 (z) := ^2u k (<T^ i 1 R<T Z2 ;n,y)e(nx). 

Theorem 5.2 ([38], Appendix B). We have ZlyZ2 \kj — Zl . Z2 for all 7 <E T, and 
the identity (13) holds. 

Remark 5.3. The quadratic form (_B,det) represents (nontrivially) if and only if 
B is split, so that U k (L; 0, y) ^ for some lattice L < B if only if B = M 2 (Q). 

5.3. Problem: computing integrals of automorphic functions. Let M be a 

positive integer and T = ro(AZ). 11 There arose in §5.1 the problem of computing 
the Petersson inner product of a fixed newform against an explicit theta series on 
r\H. This problem is a special case of a more general one to which we now turn. 

Suppose that we are given a T-invariant function F : T\M — > C. By "given," we 
mean that we know the Fourier expansion 

(14) F(tz) = ^2 a F (n, y; r)e{nx) 

of F, for each r € SL2(Z), to some large but finite precision. The main example 
to keep in mind is when F(z) = y k fi(z)f2(z) for suitably regular automorphic 
functions /1 and fi of the same weight k, at least one of which is neither holomorphic 
nor a Hecke eigenform. Other approaches are available when, say, both f\ and ji 
are holomorphic. 

We do not assume any regularity of F, but we do assume that it is absolutely 
integrable on r\H. We assume also that one can compute of (n, y;r), to a given 
precision, in essentially bounded time for large y and in time (l/y)°^ for small y. 

10 See Appendix A.l for the Fourier expansion of 6 zi ,z 2 a t other cusps. 

llr The discussion that follows applies reasonably well to any lattice V with cusps. 
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Given such an F, how can we compute the integral 

JV\M V 

In applications we might wish to compute many such integrals (say a million). 

5.4. Some candidate approaches. As motivation, we describe here some natural 
approaches to the problem stated in §5.3. 

Candidate approach 1: brute force. Suppose first that M = 1, and that the Fourier 
expansion (14) converges rapidly on the standard fundamental domain 

(15) T := {z = x + iy € H : \x\ < 1/2, \z\ > 1} 

for SLi2(Z). The latter condition holds in all applications we have in mind, so we 
henceforth refrain from stating it explicitly in this informal discussion. Then we 
can compute J F by truncating its Fourier expansion and integrating each term 
over J- . 

For general M, we can take as a fundamental domain for r\H the essentially 
disjoint union of tT over r 6 SL,2(Z)/r. 12 Then J F is the sum over r of the 
integral of F(tz) over J 7 , which we compute as above. 

Candidate approach 2: holomorphic projection. Suppose that k > 2, and that F is 
a product F(z) = y k fi(z)f2(z) where /i 6 Sk{T) is a holomorphic cusp form. Sup- 
pose also there exists 5 > such that that J2(tz) <C y~ s for all r € SL 2 (Z), y > l. 13 
Let (, ) denote the Petersson "inner product" on weight k automorphic functions on 
T, given by (hi,h 2 ) = J r ^ B y k h 1 (z)h 2 (z) whenever the integral converges ab- 
solutely. The pairing (, ) induces the structure of a finite-dimensional Hilbert space 
on iSfc(r), so there exists a unique form h 2 € Sk(X), called the holomorphic projec- 
tion of f 2 , with the property that (hi,h 2 ) — (hi,f 2 ) for all hi 6 Sfc(r). Suppose 
that the Fourier expansion of f 2 at oo reads f 2 {z) — ^2af 2 (n,y)e(nx), with the 
sum taken over n E N. Then, following Gross-Zagier [14, §IV.5], its holomorphic 
projection h 2 has the expansion h 2 (z) = ^2 ah 2 (n)e{nz) , where 

(!6) a h Jn) = — t 

' f „.k/2 r ,-2nny ~.k/2 r -2itny d x V 

Jm^ y e y e y 

By linear algebra on the finite-dimensional space <Sfc(r), we may write h 2 — c/i +h' 2 
for some complex number c and some element h 2 of 5fc(r) that is (, )-orthogonal to 
fi- The quantity J F = (/i,/ 2 ) = (/i,^) = c(/i,/i) may now be computed in a 
number of ways (see e.g. Example 5.7). 

Candidate approach 3: equidistribution of thickened horocycles. Our eventual method 
of choice (Theorem 5.6 below) is natural from several perspectives in retrospect, 
but was originally inspired by the following technique used by Holowinsky [19] in 
his work on the quantum unique ergodicity conjecture. Suppose, for simplicity, that 
M = 1. Let Y be a real parameter tending to oo. The thickened horocycle 

TZ Y :={x + iyeU: \x\ < 1/2, ye [y- 1 ^- 1 ]} 

12 In implementations, the sum over t S SL2(Z)/r is replaced by a sum over the cusps of F. 
13 This assumption holds in the setting of §5.1 provided that B M2(Q). 
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contains x Y copies of the fundamental domain T for SL2(Z)\H (see (15)), so it is 
reasonable to expect that 

If .dxdy 1 /" 2y " 1 d x y 



F* y / F(z)—^ = - / a F (0,y;l) 

Intuitively, J F should be roughly an average of the constant term cif(0, y; 1) taken 
over y sufficiently close to 0. This is made rigorous in the following lemma implicit 
in [19] (although the formulation here is our own). 

Lemma 5.4. Fix a smooth compactly- supported function h G C<?°(R^) with Mellin 
transform h A (s) — J, 2 > h(y)y~ s normalized so that h A (l) = Jgwznjg dx y 2 V ■ 
Define a norm S on the space of functions F : SL2(Z)\H — > C by S(F) — 
J zeJr y 1 / 2 \F(z)\ ^P-. Let F : SL 2 (Z)\H ->■ C be a function for which S(F) < oo, 
and let ap(0,y, 1) be its zeroth Fourier coefficient at oo as in (14). Then for each 
Y > 1, we have 

(17) / n^ = U M^M0,i,;i)^ + o(m 
jsl 2 (z)\h y r JyeR* + y \ Y 1 / 

where the implied constant depends only upon h. 

5.5. Criticism of candidate approaches. The approaches described in §5.4 
work, 1 but suffer some drawbacks, which we now describe. 

The brute force approach, which involves integrating over J-, is somewhat slow. 
In applications, we are not aware of closed form expressions for the integrals over 
T that arise, so we must compute each one afresh. A typical such integral might 
look like a linear combination of terms 

,1/2 

y'~ ~e~ e 

x 2 +y 2 >l 

over some integers b and nonnegative reals c. Caching helps, but not enough to 
carry out efficiently certain more advanced applications, such as Heegner point 
computations. 

Moreover, we have in mind the natural generalization of the problem under 
consideration to the setting of automorphic forms on GL2 over a number field, in 
which the integral over ro(M)\H with its explicit fundamental domain becomes 
replaced by an integral over a (2ri + 3r 2 ) -dimensional quotient A\(H ri x H3 2 ), or 
perhaps several copies of such quotients; here H3 is hyperbolic upper half-space. An 
explicit description of a fundamental domain for A\(H ri x H3 2 ) seems prohibitively 
complicated, even for the simplest nontrivial number fields. Even granting such a 
description, the problem of integrating an oscillating multivariate function over such 
a fundamental domain (a few million times) still seems computationally infcasible. 

The holomorphic projection approach works quite well when k is large enough, 
even over number fields. It is simpler than the other methods in the presence of 
nontrivial level because it requires knowledge of the Fourier coefficients ap{n,y;T) 
only when r = 1 (i.e., at the cusp 00). When k is small, say k — 6 or k = 4, 
the computation of the numerator of (16) becomes prohibitively slow. There is 



^We have experimented with the first two (using SAGE) for F arising as a rather general 
product of squarefree level newforms, their derivatives, and Shimizu theta series; see [25]. 
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a holomorphic projection formula for k = 2 that we have not stated here (see 
[14, §IV.6]), but it involves a delicate limiting procedure that seems to render it 
computationally ineffective. For us, the case k = 2 is important owing to its 
relevance for computing modular parametrizations of elliptic curves. 

Lemma 5.4 is convenient for certain theoretical applications, as in [19]. However, 
it is not well-suited for numerics: One must take F 1 / 2 ^> 10 R in Lemma 5.4 to 
obtain R digits of precision past the decimal point, but it may require time (not less 
than) polynomial in Y to compute the integral on the RHS of (17). For applications 
in which high precision is needed, an algorithm requiring time exponential in the 
number of digits of desired precision is infeasible. 

5.6. Our method of choice. We begin by recalling some background on Eisen- 
stein series. Let £(s) = Fr(s)£(s) be the completed Riemann zeta function, where 
r R (s) = ^- s / 2 r(s/2) and £(«) = £„ eN n~ s (Re(s) > 1). Let = {± ( 1 ™ ) : n e Z} 
be the stabilizer in SL2(Z) of the cusp oo, let 

E s (z) = J2 Im(7z) s > 1) 

7er DO \SL 2 (z) 

be the standard real-analytic Eisenstein series for SL2(Z), and let E*(z) = 2£ > {2s)E 8 {z) 
be its completion. The functions E s and E* descend to SL2(Z)\H. We collect some 
of their standard properties in the following proposition. 

Proposition 5.5 (see [20]). Let zcH. The function s H> E s (z), defined initially 
for Re(s) > 1 by an absolutely and locally uniformly convergent series, extends to 
a meromorphic function on the complex plane that is holomorphic in the half-plane 
Re(s) > 1/2 away from a simple pole at s = 1. The function s t—¥ E*{z) extends 
to a meromorphic function on the complex plane, holomorphic away from simple 
poles at s = 1 and s = 0, that satisfies the functional equation E*(z) — E\_ s {z). 

One has res s= i E s (z) — (fsL 2 (Z)XB ^p^ J an d r es s= i E*(z) — 1. The function 
E s (z) has moderate growth in both variables in the sense that for Re(s) in a fixed 
compact subset of [1/2, oo) and z = x + iy in the fundamental domain T , it satisfies 
E s (z) -C (l + lsl) ^^ 6 ^. The function s i-> E*(z) decays rapidly in vertical strips 
in the sense that if A € K, z£l and Re(s) G K are taken to lie in fixed compact 
sets, then E*(z) < (1 + |s|)- A |s(s - 1)| _1 - 

As motivation, we sketch the proof of Lemma 5.4. Let h and Y be as in its 
statement, and define E Y (z) = X) 7 er =0 \ sl 2 (z) Y~ 1 h{Y ■ lm(jz)). By unfolding, 
J E Y F gives the main term on the RHS of (17). On the other hand, Mellin inversion 
and a contour shift show that 

E Y (z) = f h\s)Y s - 1 E s (z) ^7=1 + Y-^ 2 f h A (l/2 + l t)Y* t E 1/2+tt (z) ^ 
J(2) 2tti J tm 2tt 

Since repeated partial integration implies that /i A (l/2 + it) <C (1 + |£|)~' 4 for each 
fixed A € K, we deduce from Proposition 5.5 that E Y (z) = 1 + 0(Y~ 1 / 2 y 1 / 2 ) for 
all z£ J. Integrating against F gives (17). 

We refine this argument by completing the Eisenstein series and continuing past 
the line Re(s) = 1/2, noting that unlike E s , which has infinitely many poles in the 
region < Re(s) < 1/2, the function E* is holomorphic away from s = and 8=1. 
Fix e > 0, and let H be a holomorphic function on the strip {s e C : — e < Re(s) < 
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1 + e} satisfying the conditions H(0) = 0, H(l) = 1, and H(s) < (1 + |s|)° (1) . 
For example, one can take H(s) = s. Then H(s)E* is holomorphic away from its 
simple pole at s = 1 of residue 1 and decays rapidly in the strip — e < Re(s) < 1+e, 
so for each S £ (0, e), we have the fundamental identity 



(IS) 1 = res s=1 H(s)E* s = I / - / ) H(s)E; 



'(1+5) •/(-«)/ 



2ni 

ds 



1(1+5) 27r * 

In the final step, we used the functional equation E* = E'i_ s . 

We would like to integrate (18) against F : T\M — > C and interchange the order 
of integration to obtain 

ds 



2m' 



[F = [ (H(s) - H(l - *))2£(2a) ( [ FE S 
J J (1+8) \J 

By taking S sufficiently small, we see that this interchange is justified provided that 
|F(,z)| <C y~ a for some fixed a > and all z = x + iy with y > 1. Unfolding gives 
jFE s =J y£K a (Q,y-\)y s - l d*y. 

If r is a finite index subgroup of SL^Z) and F is on r\H, then applying the 
above argument to its pushforward z h-> X)r\SL 2 (z) F( TZ ) on SL2(Z)\H yields: 

Theorem 5.6. Let V be a finite index subgroup o/SL2(Z). Le£ _F : r\H ->Cfca 
bounded measurable function satisfying F(tz) <C i/~ q for some fixed a > 0, almost 
all z = x + iy with y > 1, and a/Z r G SL 2 (Z). Lei e > 0, and let H be a holomorphic 
function on {s G C : — e < Re(s) < 1 + e} satisfying H(0) — 0, H{1) = 1, and 
H(s) < (1 + |s|)°W. Then for 5 e (0,min(a,e)), we k»e 

/ F{z) d ^ = f (H( S )-H(l-s))2Z(2s) Yl ^(0,sr) A (l- S )^. 
7r\H V 7(1+5) rer\SL 2 (z) 2 ™ 

-ffere Of(0,2/;t) is £/ie constant term of the Fourier expansion of F(tz) as in (14), 

a F (0,-;r) A (l- S ) := f a F (0, y; r)^ 1 rf x y. 

Example 5.7. We give an example in which Theorem 5.6 specializes to the approxi- 
mate functional equation applied to Rankin-Selberg L-functions. Take T = SLi2(Z), 
k = 12, and let A(z) = qf[(l - q n ) 2i = £a A (?i)q™ € S k (T) be as in §1.4. De- 
fine F : T\H -> C by the formula F(z) = y fc |A(z)| 2 . The constant term of F is 
o F (0, y; 1) = y k J2\ a A(n)\ 2 e~ ijTny , which has the Mellin transform (for Re(s) > 1) 

(n iWn , T(k-l + s) v |a A W| 2 r(fc - 1 + .s) L(A x A,s) 
n ' ' j 1 ' J (4 7 r) fc - 1 + s ^ ' (4 7 r) fe - 1 + s £(2») ' 



Thus Theorem 5.6 reads 

/|A(z)| 2 ^ = / - H(l s))2a2s)a F (0, ■; 1) A (1 - s) ^ 

SL 2 (Z)\H 2/ 7(1+5) Zm 
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where * k , H (y) = J {1+s) (H{s) - H(l - s))T( S )T(k - 1 + s)y- s For the choice 
H(s) = s, one can verify that * ktH (y) = y^-^^y/yKk^y/y) - 2JsT fc _ 1 (2^y)) 



(see Lemma A. 4). The result obtained is equivalent to what one would get by 
first expressing the Petersson norm of A in terms of the residue of L(A x A, s) 
at s = 1 via the Rankin-Selberg method, and then computing that residue via an 
approximate functional equation. The "two proofs" are essentially rearrangements 
of one another. 

Remark 5.8. Theorem 5.6 allows one to compute effectively the holomorphic pro- 
jection of certain weight k automorphic functions, even in the delicate case k = 2. 
However, it requires knowledge of Fourier expansions at all cusps. It would be in- 
teresting to have a effective method for computing holomorphic projections using 
Fourier expansions only at the cusp oo when k — 2. Conceivably one could redo 
the above argument with the non-holomorphic weight k Poincare series and its re- 
currence relation under s i— > s + 1 taking the place of the non-holomorphic weight 
Eisenstein series and its functional equation under s i-> 1 — s, but we have not 
attempted to do so. 

Remark 5.9. Theorem 5.6 is also of theoretical interest, because it gives a way to 
evaluate integrals of automorphic functions F for which one cannot satisfactorily 
control the norm S{F ) of Lemma 5.4. We exploit this property in a paper under 
preparation to bound periods of restrictions of Hilbert modular forms. 



In this section we show how the identity of Theorem 5.2 and the method for 
computing Petersson inner products afforded by Theorem 5.6 imply the identites 
stated in §3. 

A.l. Fourier expansions at various cusps. Preserve the notation of §§5.1 and 
5.2. Let F : H — > C be the function F(z) = y k f(z)8 ZlyZ2 (z), which arose as the 
integrand of (13). Owing to the identity Im(rz) = lm(z)/\cz + d\ 2 for t = ( * * d ) £ 
SLa(K), we have F(rz) = y k f\kT(z)9 ZlyZ2 \kT{z) for all r G SL 2 (R). In particular, 
Ffaz) — F(z) for all 7 e T. For each r 6 SL 2 (Z), write F(tz) =Y^ a -F(n,y; r)e{nx) 
for the Fourier expansion of F(tz), where the sum is over n in (X/w)Z with w the 
width of the cusp too. 

Suppose now that N, the level of the Eichler order R, is squarefree. Our aim in 
this subsection is to evaluate the sum, taken over r £ r\SL2(Z), of the "constant 
terms" ap(0, y;r) of F. For a lattice L < M2(K) and positive reals t and y, set 



Appendix A. Proof of formulas 



U k (L;t,y) := y k U k (L;t,y)e 



Theorem A.l. With Cf(d) as in (9), we have 



E 



a F (0,j/;r) = ^ 0/ (n) Cfidp^a^R^a^n/^y). 



rGr (d B JV)\SL 2 (Z) 



nGN d\d B N 



The proof is a calculation that occupies the remainder of this subsection. 
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Coset representatives. For each positive integer M, the natural map SL 2 (Z) — > 
SL 2 (Z/M) is surjective. By the Chinese remainder theorem, it follows that for each 
divisor d of d B N , there exists an element t c i € SL 2 (Z) such that 



Td 



1 

-1 



mod d, r d 



1 
1 



mod 



d R N 



d 



The difference between two such elements belongs to the normal subgroup 
T{d B N) = {7 G SL 2 (Z) : 7 = 1 mod d B N} of I\ The width of the cusp r d oo 
is d. For each real a; and nonzero real y, set 



n{x) 



1 a; 
1 



€ SL 2 (R), a{y) 



1/2 





-1/2 



e SL 2 (K). 



y- 

Recall that T = T (d B N) with d B N squarefree. The following lemma is well known. 

Lemma A. 2. Let d traverse the positive divisors of d B N , and let j traverse a set of 
representatives for Z/d. Then the matrices r d n{j) traverse a set of representatives 
/orr\SL 2 (Z). 

Fourier expansion of f. Let d be a divisor of d B N. The newform / £ J 7 k(d B N) 
has squarefree level. Let a/(n) be the nth Fourier coefficient of /, as in (7). By a 
well-known result of Atkin-Lehner, it transforms under Tda(d) via /|fcT<ja (d)(z) = 

^Ifid)' 1 f( z )- Since flkaid)- 1 ^) = d- k / 2 f{z/d), it follows that 

fun t\ I \ ^ ft IJ\ v(d) 
(19) f\kTd{z) = —f{z/d) = 



d - a f(d)' 



d-a f {d) 



J2a f (n)e((n/d)x)e-^ 



— 2-Kin I d)y 



Fourier expansion of zliZ2 . L e t d be a divisor of d B N. One can compute the 
Fourier expansion of Zl ,z 2 U r <i via the adelic Weil representation (see Remark B.4). 
The result obtained is that 



(20) 



'z lt Z2 \k 



;Td(z) 



fx((d, d B )) 



Ukia-iRWa* ; n/d, y)e((n/d)c 



Proof of Theorem A.l. Let d be a divisor of d B N . By (19), (20), and the definition 
(9) of Cf(d), the constant term of F(r d z) is 



d-a f (d) 



J2a f (n)U k (a- 1 R^cr Z2 ;n/d,y)e^ 



(n/d) y 



net 



Cf{d) £ a/HPk^flW^;^ y). 



d 



A direct calculation shows that ap(0, y; Tdn(j)) = a F (0,y,Td) for all integers j, so 
Lemma A. 2 gives 

E a F (0,y;T)= ^ E a F (0,y : T d n(j)) 

rero(d B JV)\SL 2 (Z) d\d B NjeZ/d 

= E«/W E c f (d)U k {a^R^o Z2] n/d,y). 

nGN d\d B N 



□ 
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A. 2. Proofs of formulas stated in §3. Let notation be as in §5.1, and set F(z) = 
y k f(z)6 Zl , Z2 ( z )- Theorem 5.2 reads 

(21) (vm) k/2 7MfB(z 2 ) = I F(z)^. 

We apply Theorem 5.6 with e = 100, H(s) = s, and S — 10, giving 

(22) / F{z) d -^L=( (2,s - l)2e(2. s ) £ a F (0, r) A (1 - .) ~. 

Suppose now that ./V is squarefree. By Theorem A.l, we then have 
(23) 

£ a F (0,-,r) A (l-s) = J2a f (n) £ c f {d)U k {a^R^a Z2 ;n/d, •)*(!-*). 

TSr\SL 2 (Z) neN dldijAr 

Lemma A. 3. For each lattice L < M%(NL) and positive real t, 

~ ds 
(2s - l)2£(2s)£/ fc (L; t, -) A (1 - s) ■— = V k (L, t). 
(10) lm 

Proof of Theorem 3.1, assuming Lemma A. 3. Combine (21), (22), (23) and Lemma 
A.3. □ 

Proof of Theorem 3.9. This is obtained from Theorem 3.1 by differentiating. We 
omit the details, but remark that the resulting formula has convincingly passed all 
numerical tests. □ 

To prove Lemma A. 3, we need the following technical lemma. 

Lemma A. 4. For each positive real x and v £ C with Rc(^) > 0, 

1 T{s)T(s + v) ds 



( s ~ I) /o\2,-u/ ^ = xK v -x(x) - K v (x). 



Proof. Mellin inversion and the integral formula / °° x s K u (2x) d x x = \F ( s ^) T ( 2 ^) 
valid for Re(s ± v) > (see [10, 6.561.16]) show that for Re(v) > -5, 

Kv{x) _i f r (^) r (^) ds _i f r( s )r( s + v) ds 



4 7(Rc M+ 5) (^/2) S 2m 2j {w) (x/2)^+- 2-ki 
For Re(i/) > -4, the identity sT(s) = T(s + 1) gives 

, . _ f t(s)t(s + u-i) d^ _ r r(s)r( s + ^) ds_ 

U - l{X) ~J m (x/2)2H-"-2 2 m - J m S (x/2)*>+» 2tT? 
Subtracting, we are done. □ 

Proof of Lemma A. 3. Since P(a) + det(a) = 2|A(a)| 2 , we see that 
U k (L;t,y) = yk ^- J2 X(a) k e-^ p ^y = yk ^ £ A( 



dot(a)=t det(a)=t 



Thus for each s £ C with Re(s) = 10, 



U k (L;t,-)*(l-s) = / U k (L;t,y)y s - L d*y 
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k r ( s + fc ) 



2 ^ v / (47r|X(a)P)-+* 

det(a)=t 



(24) = ]T e lfcar s^Wz fe (47r|X(a)|; S ), 

CtgL 

det(a)=t 

where Z k (x;s) := 2~ k ~ 1 ir s T(s + k)(x/2)~ 2s ~ k ; here the interchange of summation 
and integration is justified by absolute convergence. Since £(2s)7r s = T(s) ^2 



, m 



-2s 



it follows from Lemma A. 4 that 
(25) 

J(W) 2tti ^ J {10) (mx/2) 2s + k 2m 

= 2 x ~ k V" m k (mxK k -i(mx) - K k (mx)) 

= W k (x). 

Combining (24) and (25) (with x = 4ir\X(a)\), we conclude that 

/ (2 S -l)2£(2 S )U k (L:t,r(l- S )-^= V e* k ^ x(a )W k (4ir\X(a)\), 

J m 27TZ ^ 

det(a)=t 

which equals V k (L,t) by definition. □ 

Appendix B. Explicit Shimizu correspondence 

In this appendix, we explain how Theorem 5.2 follows from the Shimizu corre- 
spondence as pinned down by Watson [38, Thm 1]. 

B.l. Notation, measures. Let B be an indefinite rational quaternion algebra. 
Recall the notation from §§2.2 and 2.3. The discussion that follows applies to non- 
split B as well as to B = M 2 (Q). We regard B, B x , PB X , B 1 := {b G B x : 
det(fe) = 1}, and Z{B X ) = center of B x as algebraic groups over Q. 

Let A = J]' Qv = K x Q (Q = l[' p Q p ) be the adele ring of Q. For G an algebraic 

group over Q, write G v = G(Q V ). We identify G(A) = G^ x G(Q) without mention. 

Let Goo-), be the topologically connected component of G aoi which for G = 
B x or PB X is the subgroup of positive determinant elements. We fix compatible 
identifications B^ = M 2 (R) and B^ = GL 2 (R), and let = SO (2) C B^ + be 
the maximal compact connected subgroup of B^ + stabilizing iei. 

Let X = either Q„ or B v . Define (, ) : X x X -> Q„ by Q„ x 9 (x, y) n- 
and B v x B v 3 (a,/3) tr(a/? 1 ) = a/3 1 + /3a\ Let <S(X) denote the space of 
Schwarz-Bruhat functions on X. Let e = J\ e v € Hom(A/Q, S" 1 ) be the standard 
additive character, characterized by requiring 600(2;) = e(x) :— e 2 " x . Let dx be 
the Haar measure on X for which the Fourier transform T : S(X) —> S(X) given 
by J-<p(y) '■= J B (p(x)e v ((x, y}) dx satisfies TT^pix) = <p{— x). For x g X, let \x\ be 
its modulus, so that [y 1— >• xy]*dy = |arj dy. Let d x x be the measure £„(l)|x| _1 dx 
on X x , where C P (s) := (1 -p" 8 )" 1 and Coo(s) := 7r- s / 2 r(s/2). 

The measures so-defined on Q x and B x induce measures d x a on PB X and rfW 
on Bl via the short exact sequences 1 -> -> B* A (J„ x 4 1 and I -> Q* -> 
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B x — > PB X — > 1. This normalization is consistent with [38, 2.1.2]. For p ] d B 
(resp. p | ds), the measure d^a assigns volume Cp(2) — 1 (resp. (p— 1) _1 C P (2) _1 ) 
to the group of norm one units in a maximal order in B p . Let o~ z (z G H) be as in 
§3.1. Then on B^, the measure d^a is given by 

cos(6*) sin(6>) 



- dx -, — r d6 

2 \y\ 



if a 



} x-\-iy 



sin(0) cos(6>) 



The product d^a := H v dWa v of measures on B^ converges to a measure on 
B 1 (A), inducing a quotient measure on [B 1 ] := B 1 (Q)\B 1 (A) by giving B 1 (Q) the 
counting measure. We may define similarly a product measure d x a := nu^ XQ! " 
on PB X (A) inducing a quotient measure on [PB X ] := PB X (Q)\PB X (A). 



B.2. Adelization. Let R C -B M 2 (M.) be an Eichler order and real embedding 
giving rise to the lattice T = R n SL 2 (M) < SL 2 (H). Let R p = R ® z Z p and 
# = I! -Rp C S(Q). If J"(r\M) is a fundamental domain for T\M, then 

[°, ( -°!ne Ze) x «o : ze F(r\H), < < tt, k G i? x } 

is a fundamental domain for [SLa]. Fix a positive even integer A:. The following 
lemma is a consequence the strong approximation theorem. 

Lemma B.l. The identity F(a z x 1) = y k ^ 2 f(z) for z — x + iy G H induces a 
bijection / -h- F between the space of automorphic functions f of weight k onT and 
the space of smooth functions F : B^ — > C satisfying F(z^/gK oc KQ) = F{g)xk{^cxi), 
where z G Z(B X ) A , 7 G Bq , g G B^ , Koo G C B^ +) k G R x and Xk ■ #00 = 
SO(2) -> S 1 lS ^iven 6j/ ( _ c °L e e cose ) e <M - J/ /1 Fi and / 2 o F 2; tten 

1 



/ 

JIB 1 



F 1 {g)F 2 {g) dF*g 



PB> 



c(T) 



F^F^d^g 

dx dy 



y*h(z)f2{z) ., 

r\H 2/ 
-1 



wAere c(T) := 2^(2)^11^(1 "P" 1 ) IW 1 +P" 1 )- 

Remark B.2. One has c(r) = vol(r\H) (see Remark 2.3), or equivalently, volQi? 1 ]) = 
(l/2)vol([PB x ]) = 1. 

B.3. Weil representation. Let B be a quaternion algebra over Q. Let GO(B) be 
the orthogonal similitude group of the four-dimensional quadratic space (£?,det), 
and v : GO(B) -> G m the similitude factor. Let H = G(SL 2 x 0(B)) = {(g, h) G 
GL 2 x GO{B) : det(g) = u{h)}. Then H(A) acts on the space S(B(A)) = ®'S(B V ) 
of Schwarz-Bruhat functions via the Weil representation uj = ®uj v , characterized 
by (see [18, §3]) 

uj v (n(x)t(y),l)(p v {a) = \t\l /2 e v (x det(a))ip v (ta) for x G Q v ,y e Q*, 
l)Vti(a) = (~l) lv]dB F(p v (a), 

"1 



^(a) = \v(h)\ v A/4 ip v (h 1 a) for ft 



G GO(B)„. 



Here l„|d B is 1 if u divides dg and otherwise, while 



n(x) 



1 x 
1 



t(y) 
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B.4. Theta correspondence. There is a morphism p : B x x B x — > GO(B) given 
by p(6i, b 2 )(a) = b\(xb~^ , which surjects onto the connected component and satisfies 
det(6i6 2 - 1 ) = v(p(h, b 2 )). For <p G S(B A ), g G GL 2 (A) and (b x ,b 2 ) € (B x x 5 X )(A) 
with det(g) = det(6i6^" 1 ), define the theta kernel 

O v ,(g;h,b 2 ) = 2J ^(9, p( b i,b 2 ))ip(a). 

For fixed 6 X , 6 2 G B x (A) and G GL 2 (A) with det(</) = det(6i6 2 ~ 1 ), the function 
SL 2 (A) 3 g i ^ ^(53') is left-SL 2 (Q)-invariant. If F is a function on GL 2 (Q)\ GL 2 (A) 
of rapid decay modulo the center, define for b\.b 2 G -B x (A) its theta lift 

e p (F)(6 1) 6 2 )= / F{gg')e v {gg'-MM)d {1) g, 
-'gelSLa] 

where 5' G GL 2 (A) is any element for which det(g') = det(6i6 2 ~ 1 ). 

B.5. The precise lifting for newforms. Let notation be as in §5.1. Let f B ■<-» Fb 
and / o F via Lemma B.l. Define ^ = ®<^, G S{B{A)) via y> p = vol(.R x , d*^) -1 !^ 
and (^00 = it X e~^ KP , where In is the characteristic function of R p , and X and 
P are as in §3.2. This choice is consistent with [38, §2.3]. The definitions imply: 

Lemma B.3. For z,z±,z 2 G H, we have d v (a z ;a Zl ,a Z2 ) = c{T')y k / 2 9 Zl , Z2 {z). 

Remark B.4. To get the Fourier expansion for 9 ZltZ2 \Td(z) claimed in Appendix A.l, 
we use that c(Y')y k l 2 9 ZuZ2 \r d (z) = 9 v (r d a z x l;a Zl ,a Z2 ) = 9 ip (a z x t^ 1 ; <t Zi) a Z2 ) 
and apply the definition of the Weil representation at the finite places. The key 
calculation is that at a prime p for which we may identify B p = M 2 (Q p ), we have 

T\t z _ -n = p" a l, 



We omit the details, since they are discussed in [25] and [38, §2.1]. 

Theorem B.5. For b 1 ,b 2 G B x (A), we haveF B ~(b 1 )F B (b 2 )= l -^^Q ip (F)(b 1 ,b 2 ), 
where 

\\F B \\ 2 ~\l \F B \ 2 (g)d x g, \\F\\ 2 := \ f \F\ 2 (g)d x g. 

z J[PBx] z J[PGL 2 ] 

Proof. This is [38, Thm 1]. Although Watson assumes throughout his paper that N 
is squarefree, his proof of this particular result applies verbatim for general N. □ 

Proof of Theorem 5.2. By Lemma B.l, Lemma B.3 and Theorem B.5, 



(yiy2Tf B (z 1 )f B {z 2 ) = F B (a Zl )F B (a Z2 ) = Q v (F)(a Zl , a Z2 ) 



r 11 , / F(g)9 lp {g;a Zl ,a Z2 )d {1) g 
I* II -'[SL 2 ] 



I"b|| 2 

|Fb|| 2 1 f h/2 r/ \ a 1 x dxdy 

dxdy rt ||F s || 2 C (r') 

'4 J\z)o ZuZ2 \z) 

r\H 



c\ y k me ZuZ2{ z)^, c. m?c{v) 



Since f B and / are compatibly-normalized (see Definition 2.2), we have C = 1. □ 



30 paul d. nelson 

References 

[1] A. O. L. Atkin and J. Lehner. Hecke operators on T (m). Math. Ann., 185:134-160, 1970. 
[2] Henri Carayol. Sur la mauvaise reduction des courbes de Shimura. Compositio Math., 
59(2):151-230, 1986. 

[3] Henri Darmon and Robert Pollack. Efficient calculation of Stark-Heegner points via overcon- 

vergent modular symbols. Israel J. Math., 153:319-354, 2006. 
[4] Henri Darmon and Victor Rotger. Algebraic cycles and Stark-Heegner points, http:// 

www-ma2 .upc . edu/vrotger/docs/AWS2011/aws .pdf , 2011. 
[5] Lassina Dembele. Quaternionic Manin symbols, Brandt matrices, and Hilbert modular forms. 

Math. Comp., 76(258):1039-1057, 2007. 
[6] Lassina Dembele and Steve Donnelly. Computing Hilbert modular forms over fields with non- 
trivial class group. In Algorithmic number theory, volume 5011 of Lecture Notes in Comput. 

Sci., pages 371-386. Springer, Berlin, 2008. 
[7] Noam D. Elkies. Heegner point computations. In Algorithmic number theory (Ithaca, NY, 

1994), volume 877 of Lecture Notes in Comput. Sci., pages 122—133. Springer, Berlin, 1994. 
[8] Noam D. Elkies. Shimura curve computations. In Algorithmic number theory (Portland, OR, 

1998), volume 1423 of Lecture Notes in Comput. Sci., pages 1-47. Springer, Berlin, 1998. 
[9] Dorian Goldfeld. The Gauss class number problem for imaginary quadratic fields. In Heegner 

points and Rankin L-series, volume 49 of Math. Sci. Res. Inst. PubL, pages 25—36. Cambridge 

Univ. Press, Cambridge, 2004. 
[10] I. S. Gradshteyn and I. M. Ryzhik. Table of integrals, series, and products. Elsevier/ Academic 

Press, Amsterdam, seventh edition, 2007. Translated from the Russian, Translation edited 

and with a preface by Alan Jeffrey and Daniel Zwillinger, With one CD-ROM (Windows, 

Macintosh and UNIX). 

[11] Matthew Greenberg. Heegner point computations via numerical p-adic integration. In Al- 
gorithmic number theory, volume 4076 of Lecture Notes in Comput. Sci., pages 361-376. 
Springer, Berlin, 2006. 

[12] Matthew Greenberg. Heegner points and rigid analytic modular forms. ProQuest LLC, Ann 

Arbor, MI, 2006. Thesis (Ph.D.)-McGill University (Canada). 
[13] Matthew Greenberg and John Voight. Computing systems of Hecke eigenvalues associated to 

Hilbert modular forms. Math. Comp., 80(274):1071-1092, 2011. 
[14] B. Gross and D. Zagier. Heegner points and derivatives of L-series. Invent. Math., 84:225-320, 

1986. 

[15] Benedict H. Gross. Heegner points on Xa(N). In Modular forms (Durham, 1983), Ellis Hor- 
wood Ser. Math. Appl.: Statist. Oper. Res., pages 87-105. Horwood, Chichester, 1984. 

[16] Benedict H. Gross. Kolyvagin's work on modular elliptic curves. In L-functions and arith- 
metic (Durham, 1989), volume 153 of London Math. Soc. Lecture Note Ser., pages 235-256. 
Cambridge Univ. Press, Cambridge, 1991. 

[17] Paul E. Gunnells and Dan Yasaki. Hecke operators and Hilbert modular forms. In Algorithmic 
number theory, volume 5011 of Lecture Notes in Comput. Sci., pages 387-401. Springer, 
Berlin, 2008. 

[18] Michael Harris and Stephen S. Kudla. The central critical value of a triple product L-function. 

Ann. of Math. (2), 133(3):605-672, 1991. 
[19] Roman Holowinsky. Sieving for mass equidistribution. Ann. of Math. (2), 172(2):1499-1516, 

2010. 

[20] Henryk Iwaniec. Spectral methods of automorphic forms, volume 53 of Graduate Studies in 
Mathematics. American Mathematical Society, Providence, RI, second edition, 2002. 

[21] Henryk Iwaniec and Emmanuel Kowalski. Analytic number theory, volume 53 of American 
Mathematical Society Colloquium Publications. American Mathematical Society, Providence, 
RI, 2004. 

[22] Dimitar Jetchev, Kristin Lauter, and William Stein. Explicit Heegner points: Kolyva- 
gin's conjecture and non-trivial elements in the Shafarevich-Tate group. J. Number Theory, 
129(2):284-302, 2009. 

[23] David R. Kohel and Helena A. Verrill. Fundamental domains for Shimura curves. J. Theor. 

Nombres Bordeaux, 15(l):205-222, 2003. Les XXIIemes Journees Arithmetiques (Lille, 2001). 
[24] P. D. Nelson, A. Pitale, and A. Saha. Bounds for Rankin-Selberg integrals and quantum 

unique crgodicity for powerful levels. ArXiv e-prints, May 2012. 



EVALUATING MODULAR. FORMS ON SHIMURA CURVES 



31 



[25] Paul D. Nelson. Computing on Shimura curves II: Implementing the Shimizu lifting. Appendix 
C of the notes Periods and special values of L-functions, http : //swc .math. arizona. edu/aws/ 
201 l/2011PrasannaNotesProj ect . pdf , 2011. 

[26] Arnold Pizer. An algorithm for computing modular forms on Vq(N). J. Algebra, 64(2):340- 
390, 1980. 

[27] Kartik Prasanna. Computing on Shimura curves I: Expansions at CM points. Appendix B 
of the notes Periods and special values of L-functions, http://swc.math.arizona.edu/aws/ 
201 l/2011PrasannaNotesProj ect . pdf , 2011. 

[28] Hideo Shimizu. On zeta functions of quaternion algebras. Ann. of Math. (2), 81:166-193, 
1965. 

[29] Hideo Shimizu. Theta series and automorphic forms on GL2. J. Math. Soc. Japan, 24:638- 
683, 1972. 

[30] Goro Shimura. Introduction to the arithmetic theory of automorphic functions. Publications 

of the Mathematical Society of Japan, No. 11. Iwanami Shoten, Publishers, Tokyo, 1971. 

Kano Memorial Lectures, No. 1. 
[31] William Stein. Modular forms, a computational approach, volume 79 of Graduate Studies in 

Mathematics. American Mathematical Society, Providence, RI, 2007. With an appendix by 

Paul E. Gunnclls. 

[32] W. A. Stein et al. Sage Mathematics Software (Version J^.l.l). The Sage Development Team, 

2009. http://www.sagemath.org. 
[33] Marie-France Vigneras. Arithmetique des algebres de quaternions, volume 800 of Lecture 

Notes in Mathematics. Springer, Berlin, 1980. 
[34] J. Voight and J. Willis. Computing power series expansions of modular forms. ArXiv e-prints, 

April 2012. 

[35] John Voight. Computing CM points on Shimura curves arising from cocompact arithmetic 
triangle groups. In Algorithmic number theory, volume 4076 of Lecture Notes in Comput. 
Sci., pages 406-420. Springer, Berlin, 2006. 

[36] John Voight. Computing fundamental domains for Fuchsian groups. J. Theor. Nombres Bor- 
deaux, 21(2):469-491, 2009. 

[37] John Voight. Computing automorphic forms on Shimura curves over fields with arbitrary 
class number. In Algorithmic number theory, volume 6197 of Lecture Notes in Comput. Sci., 
pages 357-371. Springer, Berlin, 2010. 

[38] Thomas C. Watson. Rankin triple products and quantum chaos. arXiv . org: 0810 . 0425, 2008. 

[39] Shouwu Zhang. Heights of Heegner points on Shimura curves. Ann. of Math. (2), 153(1):27- 
147, 2001. 

EPFL, Station 8, CH-1015 Lausanne, Switzerland 
E-mail address: paul .nelsonOepf 1 . ch 



